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SUMMARY 

This paper presents a simple method for performing elastic R-curve insta- 
bility calculations. For a single material-structure combination, the calcula- 
tions can be done on some pocket calculators. On microcomputers and larger, 
it enables the development of a comprehensive program having libraries of driv- 
ing force equations for different configurations and R-curve model equations 
for different materials. The paper also presents several model equations for 
fitting to experimental R-curve data, both linear elastic and elastoplastic. 

The models are fit to data from the literature to demonstrate their viability. 


INTRODUCTION 

The R-curve is one of the most powerful concepts available to the fracture 
analyst. It is probably the best phenomenological description of the mono- 
tonic fracture process currently available. But it has not received the 
acclaim and widespread use that it deserves. Perhaps the instability calcula- 
tions are thought to be too involved or tedious. The literature contains very 
little information on instability calculations and nothing recent. Creager 
(ref. 1) presented a good graphical method. But it involves overlay transpar- 
encies, is labor-intensive, and may lack precision. In a predictive round- 
robin program (ref. 2), seven participants used the R-curve method. Only this 
author used anything more elaborate than trial -and-error to perform the insta- 
bility calculations. 

The method I used is really quite simple. For a single material-structure 
combination, the calculation can be done on some pocket calculators. On micro- 
computers and larger, a comprehensive program having libraries of driving force 
equations and R-curve model equations is possible. 

This paper will describe the method of performing the instability calcula- 
tion. It will also present several model equations for fitting to experimen- 
tal R-curve data, both linear-elastic (K,G) and elastoplastic (J). 


SYMBOLS 

a crack length 

A,B,C,D,F,H empirical coefficients in equations (3) to (7) 

G strain energy release rate 

J nonlinear crack parameter 

stress intensity factor 


K 



L.H.N.P 

W 

Y 

A 

a f 

Subscripts: 

o 

c 

R 


empirical coefficients in equations (8) to (11) 
width of specimen (or structure) 
stress intensity calibration factor 
crack extension (effective or physical, as noted) 
flow Stress 

at initial (unloaded) condition 
critical, at the instability point 
related to the material's resistance 


Specimen Notation: 

M(T) center-crack tension 

CS rectangular compact 

CLWL crack-line wedge loaded 


INSTABILITY CALCULATION METHOD 

The method was presented previously (ref. 3), but will be repeated here. 
The instability condition requires that both the magnitudes and the slopes of 
the crack driving force (G or K) curve and the material resistance curve (Gr) 
be equal at the point of instability. So we must solve two simultaneous equa- 
tions. After writing the equations in a general form and doing a little alge- 
bra, we can write the instability condition as 



G D (A ) a, + A 


n R c o i 

U = GA(A ) " 1 + 2a 
R C ' 

where 

Gpj (A) = dG R (A)/dA 

and 

a = (a/Y)(dY/da) 


For elastic R-curves A is the effective crack extension, a 0 is the ini- 
tial crack length, Y is the stress intensity calibration factor, and the sub- 
script (c) means "evaluated at the instability point." (If you prefer to work 
with the stress intensity factor K, substitute K/2K' for G/G ‘ . ) Now if we 
put in the expressions for Gr, GA, and a then, for prescribed values of a 0 
and specimen width W, A c is the K least positive root of equation (1). Any of 
several numerical methods will solve for the root, even on a programmable hand 
calculator. Once A c is known, we can calculate G c and the failure stress. 
The beauty of this equation is that the first term on the right side is a 
function of the material's R-curve only and the second is a function of the 
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structural geometry only. Equations for the second term can be obtained from 
standard handbooks and will not be discussed here. Model R-curve equations to 
be used in the first term will be discussed later. 

The effectiveness of this method is shown in appendix 11 of (ref. 2). In 
that analytical round-robin program, good predictions of fracture strengths 
were made for both simple and complex specimen geometries. 

Unfortunately this simple method is not applicable to elastic-plastic (Jr) 
calculations. As discussed in Sections (2 and 6) of (ref. 4), the problem is 
much more complicated. Material properties are an integral part of the calcu- 
lation of the driving force (J), and there appears to be no general way to sep- 
arate them from the geometrical terms. However, the model equations for Jr 
curves to be discussed later may still prove useful. 


ELASTIC R-CURVE MODELS 

To do these calculations we need an equation for the R-curve. While it 
would be nice if the equation had some physical significance, it's not abso- 
lutely necessary. Until a good theoretical analysis is available, all we need 
is a continuous equation that approximates reality, is readily differentiable, 
and fits the data. At present there is no single equation that fits all 
R-curve data. But there are several usable equations, so we have to pick the 
one that best fits the particular data. Let 1 R 1 be a generic term to indi- 
cate either Gr or Kr. 

Wang and McCabe (ref. 5) suggested fitting a polynomial to data in the 
region of the R-curve where you expect instability to occur. One could also 
describe the entire curve by a low-order (say, cubic) spline function. For 
either case we have 

R 0 = A 0 + A] A + A 2 A 2 + A 3 a3 (2) 

(note: different polynomial coefficients are used for each spline segment.) 

Broek and Vlieger (ref. 6) suggested a power curve of the form 

Rl = AA 8 (3) 

where 0<B< 1 . This particular form works very well for the more ductile materi- 
als such as 2000-series aluminums and some stainless steels. It has two unique 
features. The slope is infinite at A = 0, and there is no asymptote. In 

physical terms these imply there is no crack extension at small loads (which is 

tenable) and that a material's fracture toughness is limited only by the size 
of the specimen (which is not). So it may not be a good fundamental model, but 
it's quite good for curve fitting. 

Lower-toughness materials are usually asymptotic to a plateau value' of 
toughness. Several models fit this description. Bluhm (ref. 7) proposed an 
exponential model, 

R 2 = R 2 C1 - exp(-A/C) ] (4) 

This form has a finite slope at A = 0 and is asymptotic to R2. A somewhat 
similar model, which I will call the hyperbolic model, is 1 
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(5) 


R 3 = R*A/(D + A) 

This is actually an inverted and transposed rectangular hyperbola. Like the 

exponential model it has a finite slope at A = 0 and is asymptotic to R.j. 

Trigonometry suggests two additional asymptotic forms. They are the arctan- 
gent function 

R 4 = Rj(2/ir) arctan(irA/2F) (6) 

and the hyperbolic tangent 

R 5 = R* tanh(A/H) (7) 

The power model, equation (3), is simple enough that an illustration is 
not required. The four asymptotic models, equations (4) to (7), are shown 
schematically in figure 1(a). The first empirical coefficient (the asymptote) 
is the 'plateau value’ of fracture toughness, the second is a characteristic 
value of crack extension. Note that in all cases the initial slope is equal 
to the first empirical coefficient divided by the second. That's why the sec- 
ond coefficient was placed in the denominator of the arguments in equations 
(4), (6), and (7) rather than in the numerator. Note also that equations (4) 
to (7) are linear in the first coefficient and nonlinear in the second. Most 
nonlinear regression routines require initial estimates of the coefficients. 
These estimates can easily be made from a plot of the raw data using fig- 
ure 1(a) as a guideline. Furthermore, the fitted coefficients can easily be 
converted from one system of units to another. 

Figure (lb) shows the asymptotic models, all with the same asymptote and 
initial slope. This helps to show the inherent characteristics of each. The 
hyperbolic tangent model (curve A) has the sharpest "knee," rising quickly to 
approach the asymptote. The knee of the exponential (B) is not as sharp but 
the curve still approaches the asymptote fairly quickly. The arctangent model 
(C) rises quickly at first but then more slowly. The hyperbolic model (D) has 
the slowest approach of all. The characteristics of these models should be 
kept in mind when attempting to fit them to data. 

The low-order polynomial (or spline), equation (2), might well be the most 
accurate way to fit a single R-curve. But a comprehensive computer program 
must store the polynomial coefficients for each spline segment as well as the 
location of the knots. Interpolation for thickness or temperature effects 
could be a problem. The other R-curve models, however, only require that we 
store two constants and one equation code for each thickness and temperature. 
Interpolation for thickness or temperature should be much simpler. 


ELASTIC-PLASTIC R-CURVE MODELS 


The four asymptotic models for elastic R-curves can bg modified to repre- 
sent elastic-plastic Jr curves by replacing [R*] with [R + TA] and renaming 
the coefficients. Then equations (4) to (7) become 


'6 ' <R 6 ♦ 

T g A)[l - exp(-A/L) ] 

(8) 

7 ' (R 1 * 

T y A)A/(M + A) 

(9) 
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8 " (R 8 + 

TgA) (2/ir) arctan(ir A/2N) 

(10) 

g = <Rg + 

T g A) tanh(A/P) 

(11) 


These four models are shown schematically in figure 2. Here the first 
coefficient R* is a reference toughness value, similar to J I(; . The second 

coefficient, T, is analogous to the tearing modulus. The initial slope is 
equal to the first coefficient, R , divided by the third (L, M, N, or P). 

Using figure 2 as a guideline, all three coefficients can be estimated from a 
plot of the raw data. However, the construction is a bit more elaborate than 
for elastic R-curves. As before, the fitted coefficients may easily be con- 
verted from one system of units to another. 

At this point we have a choice. We can prescribe the initial slope as 
twice the flow stress, which is a customary assumption. Then only the parame- 
ters R* and T need to be determined empirically. Equations (8) to (11) 
then become 


6A ‘ (R 6A + T 6A A> C1 “ ex P ( - 2a fA/R 6A )] 

(8a) 

7A = <^A + T ?A A)[A/(R5 A /2o f ) + A] 

(9a) 

8A " (R 8A + Tg A AK2/ir) arctan(irafA/Rg A ) 

(10a) 

- (Rg^ + Tg^A) t & fl h ( 2 Of A / R g^ ) 

(11a) 


On first thought, prescribing the flow stress appears attractive. However, we 
still really have three parameters In the equation. We are merely fixing one 
of them. In multi parameter curve fitting, this can sometimes result in poor 
fits. Remember, too, that the flow stress is not readily measured and is only 
defined by custom. But both approaches are worthy of investigation. The coef- 
ficients of these elastic-plastic equations can also be estimated from a plot 
of the raw data (see fig. 2) and readily converted from one system of units to 
another. 


APPLICATIONS TO DATA 

To use this method and these models, we must first obtain the model equa- 
tion coefficients by fitting to experimental data. This can be done by nonlin- 
ear regression analysis. The data to follow were all fitted using the program 
MARQFIT (ref. 8) on a microcomputer. If the data used were not available in 
tabular form, published figures were enlarged xerographi cal ly and digitized. 
Fitted values of the coefficients for the data sets presented here are given 
in tables I and II. In the discussion to follow, the models of equations (8) 
to (11) will be referred to as 'modified' models and those of equations (8a) 
to (11a) as 'special' models. 


Elastic R-curves 

Figure 3 shows the power model, equation (3), fit to data for 2014-T6 alu- 
minum center-crack specimens tested at 77 K (ref. 9). Figure 4 shows the expo- 
nential model, equation (4), fit to data for 7475-T761 aluminum CLWL specimens 
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(ref. 5). Figure 5 shows the hyperbolic model, equation (5), fit to data for 
7075-T651 aluminum CLWL specimens (ref. 2). Figure 6 shows the arctangent 
model, equation (6), fit to data for 2024-T3 aluminum CLWL specimens (ref. 5). 
(Tabular data for figures 4 and 6 were obtained by private communication from 
the second author.) Figure 7 shows the hyperbolic tangent model, equation (7), 
fit to the same data as figure 4. 

These figures all show good fits, but these are not necessarily the best 
fits that could be obtained. All possible combinations of model equations and 
data sets cannot be shown for space reasons. The combinations that are pre- 
sented were selected to show that each model is a viable one. That is, there 
is at least one data set that each model can describe well. 


Elastic-Plastic R-curves 

The following figures were developed using equations (8) to (11) as 
given. That is, the initial slope (flow stress) was considered a third fit- 
ting parameter. Figure 8 shows the modified exponential model, equation (8), 
fit to data for A106C steel compact specimens at 135 °C (ref. 10). Figure 9 
shows the modified hyperbolic model, equation (9), fit to data for 2024-T351 
CLWL specimens (ref. 2). Figure 10 shows the modified arctangent model, equa- 
tion (10), fit to data for A533B steel compact specimens at 149 "C (ref. 11). 
Figure 11 shows the modified hyperbolic tangent model, equation (11), fit to 
data for A106B steel compact specimens (ref. 12). Again the figures presented 
were chosen to show that each model is a viable one. 


Elastic-Plastic R-curves (flow stress prescribed) 

Figures 12 to 15 were developed using equations (8a) to (11a). The flow 
stress was taken as the mean of the reported yield and ultimate strengths. 
Figure 12 shows the special exponential model, equation (8a), fit to the same 
data as figure 7. Figure 13 shows the special hyperbolic model fit to data 
for 5456-H117 aluminum compact specimens (ref. 13). Figure 14 shows the spe- 
cial arctangent model fit to the same data as figure 9. Figure 15 shows the 
special hyperbolic tangent model fit to the same data as figure 10. 


DISCUSSION 

All possible combinations of model equations and data sets cannot be shown 
for space reasons. The combinations that are presented were chosen to show 
that each model is a viable one. And indeed, each model is a good fit to at 
least one data set. 

For elastic R-curves, the eye is usually able to judge whether a power 
model or an asymptotic model is appropriate. But there is no obvious way to 
predict which asymptotic model will be the best fit. For the elastic-plastic 
case, there also is no way to tell in advance whether the flow stress can be 
specified or whether it should be a parameter to be fitted. However, for 4 of 
the 6 sets of Jr data, the best fit overall was obtained when the initial 
slope (flow stress) was fitted rather than prescribed. 
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CONCLUSIONS 


In summary, the instability calculation method is simple and effective. 
The R-curve model equations presented here are all viable, and at least one of 
them should fit almost any data. The empirical coefficients are easily esti- 
mated from a plot of the raw data and converted from one system of units to 
another. In combination, the method and the models make powerful tools for 
fracture analysis. 
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TABLE I. - ELASTIC R-CURVES 


Fi gure 

Equation 

Fitted equation 

3 

3 

K r = 43.92A 0 - 2175 

4 

4 

K r = 161. 3[1 - exp(-A/10.43)] 

5 

5 

K r = 50. 98A/ (0.97 18 + A) 

6 

6 

K r = 155.0(2/ir) arctan[(ir/2)(A/15.06)] 

7 

7 

K r = 158.6 tanh(A/13.76) 


TABLE II. - ELASTIC-PLASTIC R-CURVES 


Figure 

Equati on 

Fitted equation 

8 

8 

J R = (315.1 + 129.6A) [1 - exp( -A/0. 3300)] 

9 

9 

J R □ (357.7 + 3 . 600A)A/{ 5 . 1 72 + A) 

10 

10 

J R = (373.5 + 1 30 .7A) (2/ir) arctan [(ir/2) (A/0 . 3258) ] 

n 

11 

J R = (502.7 + 177. 6A) tanh(A/0.6318) 

12 

8a 

J R = (347.6 + 122.6a) [ 1 - ex p( -A/0. 4237)] 

13 

9a 

J R = (38.87 + 1 3. 00A)A/( 0.06645 + A) 

14 

10a 

J R = (396.7 + 124.2A) (2/ir) arctan (ir/2) (A/0.2493)] 

15 

11a 

J R = (553.5 + 165.2A) tanh(A/0 .8763) 
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Effective crack extension Effective crack extension 

FIGURE 1(a). - SCHEMATIC REPRESENTATION OF ELASTIC- FIGURE 1(b). - ELASTIC ASYMPTOTIC R-CURVES HAVING THE 

ASYMPTOTIC R-CURVES <C, D. F, AND H ARE COEFFI- SAME ASYMPTOTE AND INITIAL SLOPE. 

CIENTS IN EOS. (<t) TO (7)). 


Crack 

growth 

resistance 




Physical crack extension 

FIGURE 2. - SCHEMATIC REPRESENTATION OF ELASTIC- 
PLASTIC R-CURVES (L, M, N, OR P ARE COEFFICIENT? 
IN EOS. (8) TO (ID). 


Effective crack extension. A, mm. 

FIGURE 3. - POWER MODEL (EO. (3)) FIT TO DATA FOR 
2014-T6 ALUMINUM (REF. 9). 




FIGURE 4. - EXPONENTIAL MODEL (EQ. (4)) FIT TO FI6URE 5. - HYPERBOLIC MODEL (EO. (5)) FIT TO DATA 

DATA FOR 7475-T761 ALUMINUM (REF. 5). FOR 7075-T651 ALUMINUM (REF. 2). 
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Crack growth 
resistance, J„. 
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resistance. J* 
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FIGURE 6. - ARCTANGENT MODEL (EG. (6)) FIT TO 
DATA FOR 202M-T3 ALUMINUM (REF. 5). 



FIGURE 7. - HYPERBOLIC TANGENT MODEL (EQ. (7)) 
FIT TO DATA FOR 7475-T761 ALUMINUM (REF. 5). 



Effective crack extension. A , mm. 


FIGURE 8. - MODIFIED EXPONENTIAL MODEL (EQ. (8)) 
FIT TO DATA FOR A106C STEEL (REF. 10). 


Crack growth 
resistance. 
KJ/m 2 



Physical crack extension. A . mm 

FIGURE 9. - MODIFIED HYPERBOLIC MODEL (EQ. (9)) 
FIT TO DATA FOR 2024-T351 ALUMINUM (REF. 2). 



Physical crack extension, A . mm. 

FIGURE 10. - MODIFIED ARCTANGENT MODEL 
(EQ. (10)) FIT TO DATA FOR A535B STEEL 
(REF. 11). 



FIGURE 11. - MODIFIED HYPERBOLIC TANGENT MODEL 
(EQ. (ID) FIT TO DATA FOR A106B STEEL 
(REF. 12). 
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Physical crack extention. A . mm. 


FIGURE 12. - SPECIAL EXPONENTIAL MODEL <EQ. (8a)) 
FIT TO DATA FOR A106C STEEL (REF. 10). 


FIGURE 13. - SPECIAL HYPERBOLIC MODEL (EQ. (9a)) 
FIT TO DATA FOR 5456-H117 ALUMINUM (REF. 13). 



Physical crack extension. A . rrm 

FIGURE 1M. - SPECIAL ARCTANGENT MODEL (EO. (10a)) 
FIT TO DATA FOR A533B STEEL (REF. 11). 



Physical crack extension A . mm. 


FIGURE 15. - SPECIAL HYPERBOLIC TANGENT MODEL 
(EQ. (11 a)) FIT TO DATA FOR A533B STEEL 
(REF. 11). 
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